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Abstract 

Uncertainties on low-x PDFs are crucial for the standard model benchmark 
processes of W and Z production at the LHC. The current level of PDF un- 
certainty is critically reviewed and the possibility of reducing this uncertainty 
using early LHC data is investigated taking into account realistic expectations 
for measurement accuracy, kinematic cuts and backgrounds. 

1 Introduction 

At leading order (LO), W and Z production occur by the process, qq — ► W/Z, and the momentum 
fractions of the partons participating in this subprocess are given by, x\p = ^ ex P{^y)> where M is 
the centre of mass energy of the subprocess, M = My/ or Mz, y/s is the centre of mass energy of 
the reaction (i/s = 14 TeV at the LHC) and y = \}n ^^l gives the partem rapidity. The kinematic 
plane for LHC parton kinematics is shown in Fig. ^ Thus, at central rapidity, the participating partons 
have small momentum fractions, x ~ 0.005. Moving away from central rapidity sends one parton to 
lower x and one to higher x, but over the measurable rapidity range, \y\ < 2.5, x values remain in 
the range, 10~ 4 < x < 0.1. Thus, in contrast to the situation at the Tevatron, valence quarks are not 
involved, the scattering is happening between sea quarks. Furthermore, the high scale of the process 
Q 2 = M 2 ~ 10, 000 GeV 2 ensures that the gluon is the dominant parton, see Fig. Q so that these sea 
quarks have mostly been generated by the flavour blind g — > qq splitting process. Thus the precision of 
our knowledge of W and Z cross-sections at the LHC is crucially dependent on the uncertainty on the 
momentum distribution of the gluon. 

HERA data have dramatically improved our knowledge of the gluon, as illustrated in Fig. 12 which 
shows W and Z rapidity spectra predicted from a global PDF fit which does not include the HERA data, 
compared to a fit including HERA data. The latter fit is the ZEUS-S global fit [1], whereas the former is a 
fit using the same fitting analysis but leaving out the ZEUS data. The full PDF uncertainties for both fits 
are calculated from the error PDF sets of the ZEUS-S analysis using LHAPDF [2] (see the contribution 
of M.Whalley to these proceedings). The predictions for the W/Z cross-sections, decaying to the lepton 
decay mode, are summarised in Table ^ The uncertainties in the predictions for these cross-sections 
have decreased from ~ 16% pre-HERA to ~ 3.5% post-HERA. The reason for this can be seen clearly 
in Fig. |3] where the sea and gluon distributions for the pre- and post-HERA fits are shown for several 
different Q 2 bins, together with their uncertainty bands. It is the dramatically increased precision in the 
low-x gluon PDF, feeding into increased precision in the low-x sea quarks, which has led to the increased 
precision on the predictions for W/Z production at the LHC. 

Further evidence for the conclusion that the uncertainties on the gluon PDF at the input scale 
(Qq = 7 GeV 2 , for ZEUS-S) are the major contributors to the uncertainty on the W/Z cross-sections at 
Q 2 = M\y(Mz), comes from decomposing the predictions down into their contributing eigenvectors. 
Fig|4]shows the dominant contributions to the total uncertainty from eigenvectors 3, 7, and 11 which are 
eigenvectors which are dominated by the parameters which control the low-x, mid-x and high-x, gluon 
respectively. 

The post-HERA level of precision illustrated in Fig.|2]is taken for granted in modern analyses, such 
that W/Z production have been suggested as 'standard-candle' processes for luminosity measurement. 
However, when considering the PDF uncertainties on the Standard Model (SM) predictions it is necessary 
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Fig. 2: LHC W + ,W~ , Z rapidity distributions and their PDF uncertainties (the full line shows the central value and the dashed 
lines show the spread of the uncertainty): Top Row: from the ZEUS-S global PDF analysis not including HERA data; left plot 
W + ; middle plot VK"; right plot Z: Bottom Row: from the ZEUS-S global PDF analysis including HERA data; left plot W + ; 
middle plot W~ ; right plot Z 




Fig. 3: Sea (xS) and gluon (xg) PDFs at various Q 2 : left plot; from the ZEUS-S global PDF analysis not including HERA 
data; right plot: from the ZEUS-S global PDF analysis including HERA data. The inner cross-hatched error bands show 
the statistical and uncorrelated systematic uncertainty, the outer error bands show the total uncertainty including experimental 
correlated systematic uncertainties, normalisations and model uncertainty. 




Fig. 4: LHC W rapidity distributions and their PDF uncertainties due to the eigenvectors 3,7 and 1 1 of the ZEUS-S analysis. 



PDF Set a{W+).B(W+ -> l+vj) a(W-).B(W- -> l~vi) a(Z).B(Z ^ l + l~) 

ZEUS-SnoHERA 10.63 ± 1.73 nb 7.80 ± 1.18 nb 1.69 ± 0.23 nb 

ZEUS-S 12.07 ± 0.41 nb 8.76 ± 0.30 nb 1.89 ± 0.06 nb 

CTEQ6.1 11.66 ± 0.56 nb 8.58 ± 0.43 nb 1.92 ± 0.08 nb 

MRST01 11.72 ± 0.23 nb 8.72 ± 0.16 nb 1.96 ± 0.03 nb 




Fig. 5: LHC W + rapidity distributions and their PDF uncertainties: left plot, ZEUS-S PDFs; middle plot, CTEQ6.1 PDFs; 
right plot: MRST01 PDFs. 



not only to consider the uncertainties of a particular PDF analysis, but also to compare PDF analyses. 
Fig.|5]compares the predictions for W + production for the ZEUS-S PDFs with those of the CTEQ6.1 [3] 
PDFs and the MRST01 [4] PDFs 1 . The corresponding W + cross-sections, for decay to leptonic mode 
are given in Tabled Comparing the uncertainty at central rapidity, rather than the total cross-section, we 
see that the uncertainty estimates are rather larger: 5.2% for ZEUS-S; 8.7% for CTEQ6.1M and about 
3.6% for MRST01. The difference in the central value between ZEUS-S and CTEQ6.1 is 3.5%. Thus 
the spread in the predictions of the different PDF sets is comparable to the uncertainty estimated by the 
individual analyses. Taking all of these analyses together the uncertainty at central rapidity is about 8%. 

Since the PDF uncertainty feeding into the W + , W~ and Z production is mostly coming from 
the gluon PDF, for all three processes, there is a strong correlation in their uncertainties, which can be 
removed by taking ratios. Fig.[6]shows the W asymmetry 

A w = (W + - W~)/(W + + W~). 

for CTEQ6.1 PDFs, which have the largest uncertainties of published PDF sets. The PDF uncertainties on 
the asymmetry are very small in the measurable rapidity range. An eigenvector decomposition indicates 
that sensitivity to high-x it and d quark flavour distributions is now evident at large y. Even this residual 
flavour sensitivity can be removed by taking the ratio 

A zw = Z/{W + + W~) 

as also shown in Fig.|6] This quantity is almost independent of PDF uncertainties. These quantities have 
been suggested as benchmarks for our understanding of Standard Model Physics at the LHC. However, 
whereas the Z rapidity distribution can be fully reconstructed from its decay leptons, this is not possible 
for the W rapidity distribution, because the leptonic decay channels which we use to identify the W's 

'MRSTOl PDFs are used because the full error analysis is available only for this PDF set. 



Fig. 6: Predictions for W, Z production at the LHC from the CTEQ6. 1 PDFs. Top row: left plot, the W asymmetry, Aw', right 
plot, the ratio, Azw- Bottom row: left plot, decay e + rapidity spectrum; middle plot, decay e~ rapidity spectrum; right plot, 
lepton asymmetry, A e 

have missing neutrinos. Thus we actually measure the W's decay lepton rapidity spectra rather than the 
W rapidity spectra. The lower half of Fig. [5] shows the rapidity spectra for positive and negative leptons 
from W + and W~ decay and the lepton asymmetry, 

A l = (i+-r)/{i+ + r). 

A cut of, pu > 25 GeV, has been applied on the decay lepton, since it will not be possible to trigger 
on leptons with small p^. A particular lepton rapidity can be fed from a range of W rapidities so that 
the contributions of partons at different x values is smeared out in the lepton spectra, but the broad 
features of the W spectra and the sensitivity to the gluon parameters remain. The lepton asymmetry 
shows the change of sign at large y which is characteristic of the V — A structure of the lepton decay. 
The cancellation of the uncertainties due to the gluon PDF is not so perfect in the lepton asymmetry as 
in the W asymmetry. Nevertheless in the measurable rapidity range sensitivity to PDF parameters is 
small. Correspondingly, the PDF uncertainties are also small ( 4%) and this quantity provides a suitable 
Standard Model benchmark. 

In summary, these preliminary investigations indicate that PDF uncertainties on predictions for the 
W, Z rapidity spectra, using standard PDF sets which describe all modern data, have reached a precision 
of ~ 8%. This may be good enough to consider using these processes as luminosity monitors. The 
predicted precision on ratios such as the lepton ratio, A\, is better (~ 4%) and this measurement may 
be used as a SM benchmark. It is likely that this current level of uncertainty will have improved before 
the LHC turns on- see the contribution of C. Gwenlan to these proceedings. The remainder of this 
contribution will be concerned with the question: how accurately can we measure these quantities and 
can we use the early LHC data to improve on the current level of uncertainty. 

2 k-factor and PDF re-weighting 

To investigate how well we can really measure W production we need to generate samples of Monte- 
Carlo (MC) data and pass them through a simulation of a detector. Various technical problems arise. 
Firstly, many physics studies are done with the HERWIG (6.505) [5], which generates events at LO with 



Fig. 7: Top Row: W rapidity and pt spectra for events generated with HERWIG + k-Factors (full line), compared to those 
generated by MC@NLO (dashed line); left plot W rapidity; middle plot W~ rapidity; right plot W~ pt- Bottom row: 
the fractional differences of the spectra generated by HERWIG + k-factors and those generated by MC@NLO. The full line 
represents the weighted mean of these difference spectra and the dashed lines show its uncertainty 

parton showers to account for higher order effects. Distributions can be corrected from LO to NLO by 
k-factors which are applied as a function of the variable of interest. The use of HERWIG is gradually 
being superceded by MC@NLO (2.3) [6] but this is not yet implemented for all physics processes. Thus 
it is necessary to investigate how much bias is introduced by using HERWIG with k-factors. Secondly, to 
simulate the spread of current PDF uncertainties, it is necessary to run the MC with all of the eigenvector 
error sets of the PDF of interest. This would be unreasonably time-consuming. Thus the technique of 
PDF reweighting has been investigated. 

One million W — > eu e events were generated using HERWIG (6.505). This corresponds to 43 
hours of LHC running at low luminosity, 10/6 -1 . These events are split into W + and W~ events 
according to their Standard Model cross-section rates, 58%: 42% (the exact split depends on the input 
PDFs). These events are then weighted with k-factors, which are analytically calculated as the ratio 
of the NLO to LO cross-section as a function of rapidity for the same input PDF [7]. The resultant 
rapidity spectra for W + , W~ are compared to rapidity spectra for ~ 107, 700 events generated using 
MC@NLO(2.3) in FigEf. The MRST02 PDFs were used for this investigation. The accuracy of this 
study is limited by the statistics of the MC@NLO generation. Nevertheless it is clear that HERWIG with 
k-factors does a good job of mimicking the NLO rapidity spectra. However, the normalisation is too high 
by 3.5%. This is not suprising since, unlike the analytic code, HERWIG is not a purely LO calculation, 
parton showering is also included. This normalisation difference is not too crucial since in an analysis 
on real data the MC will only be used to correct data from the detector level to the generator level. For 
this purpose, it is essential to model the shape of spectra to understand the effect of experimental cuts 
and smearing but not essential to model the overall normalisation perfectly. However, one should note 
that HERWIG with k-factors is not so successful in modelling the shape of the p t spectra, as shown in 
the right hand plot of Fig. This is hardly surprising, since at LO the W have no p t and non-zero 
Pt for HERWIG is generated by parton showering, whereas for MC@NLO non-zero p t originates from 
additional higher order processes which cannot be scaled from LO, where they are not present. 

Suppose we generate W events with a particular PDF set: PDF set 1. Any one event has the 
hard scale, Q 2 = M\y, and two primary partons of flavours flav\ and flav2, with momentum fractions 
x\, X2 according to the distributions of PDF set 1. These momentum fractions are applicable to the hard 
process before the parton showers are implemented in backward evolution in the MC. One can then 
evaluate the probability of picking up the same flavoured partons with the same momentum fractions 

2 In MC@NLO the hard emissions are treated by NLO computations, whereas soft/collinear emissions are handled by the 
MC simulation. In the matching procedure a fraction of events with negative weights is generated to avoid double counting. 
The event weights must be applied to the generated number of events before the effective number of events can be converted to 
an equivalent luminosity. The figure given is the effective number of events. 



Fig. 8: Left side: W~ (left) and W (right) rapidity spectra, for events generated with MRST02 PDFs reweighted to CTEQ6. 1 
PDFs (full line), compared to events generated directly with CTEQ6.1 PDFs (dashed line). The fractional difference between 
these spectra are also shown beneath the plots. The full line represents the weighted mean of these difference spectra and the 
dashed lines show its uncertainty. Right side: the same for p t spectra. 

from an alternative PDF set, PDF set 2, at the same hard scale. Then the event weight is given by 

urn? • fpDF 2 (xi, flavi, Q 2 ).fppF 2 (x 2 , flav 2 , Q 2 ) 

PDFre - weight = (1) 

fpDFi (xi , flavi , Q ) -iPDFi (X2 , nav 2 , Q z ) 

where xfpDF{x,flav,Q 2 ) is the parton momentum distribution for flavour, flav, at scale, Q 2 , and 
momentum fraction, x. Fig. [8] compares the W + and W~ spectra for a million events generated using 
MRST02 as PDF set 1 and re- weighting to CTEQ6.1 as PDF set 2, with a million events which are di- 
rectly generated with CTEQ6. 1. Beneath the spectra the fractional difference between these distributions 
is shown. These difference spectra show that the reweighting is good to better than 1%, and there is no 
evidence of a y dependent bias. This has been checked for reweighting between MRST02, CTEQ6.1 and 
ZEUS-S PDFs. Since the uncertainties of any one analysis are similar in size to the differences between 
the analyses it is clear that the technique can be used to produce spectra for the eigenvector error PDF 
sets of each analysis and thus to simulate the full PDF uncertainties from a single set of MC generated 
events. Fig. [8] also shows a similar comparison for p t spectra. 

3 Background Studies 

To investigate the accuracy with which W events can be measured at the LHC it is necessary to make 
an estimate of the importance of background processes. We focus on W events which are identified 
through their decay to the W — » e u e channel. There are several processes which can be misidentified 
as W — > eu e . These are: W — ► rv T , with r decaying to the electron channel; Z — > t + t~ with at least 
one r decaying to the electron channel (including the case when both r's decay to the electron channel, 
but one electron is not identified); Z — ► e + e~ with one electron not identified. We have generated one 
million events for each of these background processes, using HERWIG and CTEQ5L, and compared 
them to one million signal events generated with CTEQ6.1. We apply event selection criteria designed 
to eliminate the background preferentially. These criteria are: 

• ATLFAST cuts (see Sec.|5} 

• pseudorapidity, \rj\ < 2.4, to avoid bias at the edge of the measurable rapidity range 

• pte > 25 GeV, high p t is necessary for electron triggering 

• missing E t > 25 GeV, the v e in a signal event will have a correspondingly large missing E t 

• no reconstructed jets in the event with p t > 30 GeV, to discriminate against QCD background 



Table 2: Reduction of signal and background due to cuts 



Cut 


W -> eu e 


Z -> T+T- 

e + e~ 


Z -> e+e' 

e + e~ 


W -> TV T 


ATLFAST cuts 

M < 2.4 
pte > 25 GeV 
Ptmtss > 25 GeV 
No jets withP t > 30 GeV 

p re COl l < 2 Ge y 


382,902 264,415 
367,815 255,514 
252,410 194,562 
212,967 166,793 
187,634 147,415 
159,873 125,003 


5.5% 7.9% 
5.5% 7.8% 
0.6% 0.7% 
0.2% 0.2% 
0.1% 0.1% 
0.1% 0.1% 


34.7% 50.3% 
34.3% 49.4% 
12.7% 16.2% 
0.1% 0.2% 
0.1% 0.1% 
0.0% 0.0% 


14.8% 14.9% 
14.7% 14.8% 
2.2% 2.3% 
1.6% 1.6% 
1.2% 1.2% 
1.2% 1.2% 


• recoil on the 


transverse plane p r t ecml < 20 GeV, to discriminate against QCD background 



Table El gives the percentage of background with respect to signal, calculated using the known relative 
cross-sections of these processes, as each of these cuts is applied. After, the cuts have been applied the 
background from these processes is negligible. However, there are limitations on this study from the fact 
that in real data there will be further QCD backgrounds from 2 — > 2 processes involving q, q, g in which 
a final state -kq — > 77 decay mimics a single electron. A preliminary study applying the selection criteria 
to MC generated QCD events suggests that this background is negligible, but the exact level of QCD 
background cannot be accurately estimated without passing a very large number of events though a full 
detector simulation, which is beyond the scope of the current contribution. 

4 Charge misidentification 

Clearly charge misidentification could distort the lepton rapidity spectra and dilute the asymmetry A\. 

A — F + + F~ 
Atrue ~ 1 _ p- _ p+ 

where A raw is the measured asymmetry, A true is the true asymmetry, F~ is the rate of true misiden- 
tified as e + and F + is the rate of true e + misidentified as e~. To make an estimate of the importance of 
charge misidentification we use a sample of Z — > e + e~ events generated by HERWIG with CTEQ5L 
and passed through a full simulation of the ATLAS detector. Events with two or more charged electro- 
magnetic objects in the EM calorimeter are then selected and subject to the cuts; \r)\ < 2.5, pt e > 25 
GeV, as usual and, E/p < 2, for bremsstrahlung rejection. We then look for the charged electromagnetic 
pair with invariant mass closest to Mz and impose the cut, 60 < Mz < 120 GeV. Then we tag the 
charge of the better reconstructed lepton of the pair and check to see if the charge of the second lepton is 
the same as the first. Assuming that the pair really came from the decay of the Z this gives us a measure 
of charge misidentification. Fig [9] show the misidentification rates F + , F~ as functions of pseudorapid- 
ity 3 . These rates are very small. The quantity Ai, can be corrected for charge misidentification applying 
Barlow's method for combining asymmetric errors [8]. The level of correction is 0.3% in the central 
region and 0.5% in the more forward regions. 

5 Compare events at the generator level to events at the detector level 

We have simulated one million signal, W — ► ev e , events for each of the PDF sets CTEQ6.1, MRST2001 
and ZEUS-S using HERWIG (6.505). For each of these PDF sets the eigenvector error PDF sets have 
been simulated by PDF reweighting and k-factors have been applied to approximate an NLO generation. 
The top part of Fig. ^0] shows the and A\ spectra at this generator level, for all of the PDF sets 
superimposed. The events are then passed through the ATLFAST fast simulation of the ATLAS detector. 
This applies loose kinematic cuts: \rj\ < 2.5, pt e > 5 GeV, and electron isolation criteria. It also smears 

3 These have been corrected for the small possibility that the better reconstructed lepton has had its charge misidentified as 
follows. In the central region, \r]\ < 1, assume the same probability of misidentification of the first and second leptons, in the 
more forward regions assume the same rate of first lepton misidentification as in the central region. 




Fig. 9: The rates of charge misidenttification as a function of rapidity for e misidentified as e + (left), e + misidentifed as e 
(fight). 
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Fig. 10: Top row: e~ , e + and A e rapidity spectra for the lepton from the W decay, generated using HERWIG + k factors and 
CTEQ6.1 (red), ZEUS-S (green) and MRST2001 (black) PDF sets with full uncertainties. Bottom row: the same spectra after 
passing through the ATLFAST detector simulation and selection cuts. 



the 4-momenta of the leptons to mimic momentum dependent detector resolution. We then apply the 
selection cuts described in Sec.|3] The lower half of Fig.Q^J shows the e^- and Ai spectra at the detector 
level after application of these cuts, for all of the PDF sets superimposed. The level of precision of each 
PDF set, seen in the analytic calculations of Fig. |5j is only slightly degraded at detector level, so that 
a net level of PDF uncertainty at central rapidity of ~ 8% is maintained. The anticipated cancellation 
of PDF uncertainties in the asymmetry spectrum is also observed, within each PDF set, and the spread 
between PDF sets suggests that measurements which are accurate to better than ~ 5% could discriminate 
between PDF sets. 

6 Using LHC data to improve precision on PDFs 

The high cross-sections for W production at the LHC ensure that it will be the experimental systematic 
errors, rather than the statistical errors, which are determining. We have imposed a random 4% scat- 
ter on our samples of one million W events, generated using different PDFs, in order to investigate if 
measurements at this level of precision will improve PDF uncertainties at central rapidity significantly 
if they are input to a global PDF fit. Fig. shows the e + and rapidity spectra for events generated 
from the ZEUS-S PDFs (\r)\ < 2.4) compared to the analytic predictions for these same ZEUS-S PDFs. 
The lower half of this figure illustrates the result if these events are then included in the ZEUS-S PDF 
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Fig. 11: Top row: e + and e rapidity spectra generated from ZEUS-S PDFs compared to the analytic prediction using ZEUS- 
S PDFs. Bottom row: the same lepton rapidity spectra compared to the analytic prediction AFTER including these lepton 
pseudo-data in the ZEUS-S PDF fit. 



fit. The size of the PDF uncertainties, at y = 0, decreases from 5.8% to 4.5%. The largest improvement 
is in the PDF parameter X g controlling the low-x gluon at the input scale, Q^: xg(x) ~ x Xs > at low-x, 
X g = -0.199 ± 0.046, before the input of the LHC pseudo-data, compared to, X g = -0.196 ± 0.029, 
after input. Note that whereas the relative normalisations of the e + and e~ spectra are set by the PDFs, 
the absolute normalisation of the data is free in the fit so that no assumptions are made on our ability to 
measure luminosity. Secondly, we repeat this procedure for events generated using the CTEQ6.1 PDFs. 
As shown in Fig. El the cross-section for these events is on the lower edge of the uncertainty band of 
the ZEUS-S predictions. If these events are input to the fit the central value shifts and the uncertainty de- 
creases. The value of the parameter X g becomes, X g = —0.189 ± 0.029, after input of these pseudo-data. 
Finally to simulate the situation which really faces experimentalists we generate events with CTEQ6.1, 
and pass them through the ATLFAST detector simulation and cuts. We then correct back from detector 
level to generator level using a different PDF set- in this cases the ZEUS-S PDFs- since in practice we will 
not know the true PDFs. Fig. [T3l shows that the resulting corrected data look pleasingly like CTEQ6.1, 
but they are more smeared. When these data are input to the PDF fit the central values shift and errors 
decrease just as for the perfect CTEQ6.1 pseudo-data. The value of X g becomes, A = —0.181 ± 0.030, 
after input of these pseudo-data. Thus we see that the bias introduced by the correction procedure from 
detector to generator level is small compared to the PDF uncertainty. 



7 Conclusions and a warning: problems with the theoretical predictions at small-x? 

We have investigated the PDF uncertainty on the predictions for W and Z production at the LHC, taking 
into account realistic expectations for measurement accuracy and the cuts on data which will be needed 
to identify signal events from background processes. We conclude that at the present level of PDF 
uncertainty the decay lepton asymmetry, Ai, will be a useful standard model benchmark measurement, 
and that the decay lepton spectra can be used as a luminosity monitor which will be good to ~ 8%. 
However, we have also investigated the measurement accuracy necessary for early measurements of 
these decay lepton spectra to be useful in further constraining the PDFs. A systematic measurement 





Fig. 12: Top row: e + and e rapidity spectra generated from CTEQ6. 1 PDFs compared to the analytic prediction using ZEUS- 
S PDFs. Bottom row: the same lepton rapidity spectra compared to the analytic prediction AFTER including these lepton 
pseudo-data in the ZEUS-S PDF fit. 
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Fig. 13: Top row: e + and rapidity spectra generated from CTEQ6. 1 PDFs, which have been passed through the ATLFAST 
detector simulation and corrected back to generator level using ZEUS-S PDFs, compared to the analytic prediction using 
ZEUS-S PDFs. Bottom row: the same lepton rapidity spectra compared to the analytic prediction AFTER including these 
lepton pseudo-data in the ZEUS-S PDF fit. 




Fig. 14: LHC W + , W~, Z rapidity distributions for the MRST03 PDFs: left plot W + ; middle plot W~; right plot Z 

error of less than ~ 4% would provide useful extra constraints. 

However, a caveat is that the current study has been performed using standard PDF sets which are 
extracted using NLO QCD in the DGLAP [9] formalism. The extension to NNLO is straightforward, 
giving small corrections ~ 1%. PDF analyses at NNLO including full accounting of the PDF uncertain- 
ties are not extensively available yet, so this small correction is not pursued here. However, there may 
be much larger uncertainties in the theoretical calculations because the kinematic region involves low-x. 
There may be a need to account for ln(l/x) resummation (first considered in the BFKL [10] formalism) 
or high gluon density effects. See reference [1 1] for a review. 

The MRST group recently produced a PDF set, MRST03, which does not include any data for 
x < 5 x 10~ 3 . The motivation behind this was as follows. In a global DGLAP fit to many data sets there 
is always a certain amount of tension between data sets. This may derive from the use of an inappropriate 
theoretical formalism for the kinematic range of some of the data. Investigating the effect of kinematic 
cuts on the data, MRST found that a cut, x > 5 x 1CT 3 , considerably reduced tension between the 
remaining data sets. An explanation may be the inappropriate use of the DGLAP formalism at small-x. 
The MRST03 PDF set is thus free of this bias BUT it is also only valid to use it for x > 5 x 1CT 3 . 
What is needed is an alternative theoretical formalism for smaller x. However, the MRST03 PDF set 
may be used as a toy PDF set, to illustrate the effect of using very different PDF sets on our predictions. 
A comparison of Fig. El with Fig. |2]or Fig.|5]shows how different the analytic predictions are from the 
conventional ones, and thus illustrates where we might expect to see differences due to the need for an 
alternative formalism at small-x. 
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